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ABSTRACT 

The MAGIC collaboration observed BL Lacertae for 22.2 hr during 2005 August to December and for 26 hr 
during 2006 July to September. The source is the historical prototype and eponym of a class of low-frequency 
peaked BL Lacertae (LBL) objects. A very high energy (VHE) 7-ray signal was discovered with a 5.1 a excess 
in the 2005 data. Above 200 GeV, an integral flux of (0.6 ± 0.2) x 10 -11 cm -2 s -1 was measured, corresponding 
to approximately 3% of the Crab flux. The differential spectrum between 150 and 900 GeV is rather steep with 
a photon index of -3.6 ± 0.5. The light curve shows no significant variability during the observations in 2005. 
For the first time a clear detection of VHE 7-ray emission from an LBL object was obtained with a signal below 
previous upper limits. The 2006 data show no significant excess. This drop in flux follows the observed trend 
in optical activity. 

Subject headings: BL Lac objects: individual (BL Lacertae) - gamma rays: observations 



1. INTRODUCTION 

BL Lacertae (1ES2200+420, z=0.069. lMiiler et all 1978b is 
the historical prototype of a class of powerful 7-ray emitters: 
"BL Lac objects". It belongs to an active galactic nucleus 
(AGN) subclass in which the jet is aligned very close to our 
line of sight. The mass of the supermassive black hole in the 
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center is estimated to be ~ 10 8 M© (IWu et al.ll2QQ2h . Electro- 
magnetic emission from this class of sources can be observed 
from the radio up to very high energy (VHE) 7-ray s (E > 100 
GeV), with spectral energy distributions (SEDs) characterized 
by a two-bump structure. So far, most of the measured SEDs 
could be interpreted using scenarios of a simple leptonic ori- 
gin. The lower energy bump is produced by synchrotron ra- 
diation of relativistic electrons, while the higher energy bump 
originates from inverse Compton (IC) scattering of the same 
electron population accelerated in the jet. The target pho- 
tons for IC scattering could be either synchrotron photons 
(synchrotron self-Compton scattering, SSC) or external pho- 
tons either from the broad-line emission region or from the 
accretion disc (external inverse Compton scattering, EC). If 
both photon populations contribute to the IC emission, the 
corresponding IC component may h ave a double local-peak 
structure (e.g. Ghisellin i et al.|[l998l) . When the synchrotron 
emission peak is located in the sub-millimeter to optical 
band, the objects are classified as "Low-frequency peaked BL 
Lacs" (LBLs), whereas in "High-frequency peaked BL Lacs" 
(HB Ls) the synchrotron peak is located at UV to X -ray ener- 
gies dPadovani & Giommi|[l995l: iFossati et al.lll998h . BL Lac 
objects often show a strong flux variability down to timescales 
of a few minutes for the h ighest energy range dGaidos et al.l 
119961: lAlbert et al.l l2007d) . It should be noted that other 
models, based e.g. on the accele ration of hadrons (Mannheim 
l!993l:lMucke & Protherod2001h . could also explain the SEDs 
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Fig. 1 . — Sky map of the region around the position of BL Lacertae (black 
cross) for reconstructed 7 events with > 350 photo-electrons (corresponding 
to an energy threshold of about 200 GeV) from the 2005 observations. 



of BL Lac objects. 

Up to now, there are 13 well-established extra-galactic VHE 
7-ray sources. Twelve are classified as HBL and one is the 
Fanaroff-Riley class I radio galaxy, M87. BL Lacertae is clas- 
sified as an LB L object with a synchr otron peak frequency of 
2.2 x 10 14 Hz (ISambruna et al.ll 19991) . and is one of the best- 
studied objects in the various energy bands. In 1998 and 2001, 
VHE 7-ray de tection of the LBLs 3C66A and BL Lacertae 
were claimed (Nes hpor et al.lll998ll200ll) . Subsequent obser- 
vations with more sensitive instruments have not, however, re- 
vealed any evidence for VHE 7-ray emission consistent with 
these claims. No VHE 7-ray emission from any LBLs has 
p reviously bee n confirmed. 

IVillata et all (12004)) have presented long-term light curves 
in optical and radio emission over 30 years and reported cross- 
correlation between the optical light curve and radio hardness 
ratio with some delay in radio emission. They have claimed 
evidence of an ~8 yr periodicity in radio but less evidence 
in optical. Severa l authors (e.g. lHagen-Thorn et al.1 1 19971: 
Stirli ng etail 120031) have also reported periodic and quasi- 
periodic variations in the optical and radio light curves. 

Gamma-ray observations by EGRET resulted only in an up- 
per limit of 1.4 x 10" 7 cm" 2 s" 1 until 1995. Later, EGRET ob- 



served 7-ray s above 100 MeV at a flu x level of (4.0±1.2) x 
10~ 7 cm -2 s -1 with 4.4 a significance dCatanese et al.lll997|) . 
During an optical outburst in 1997, a 7-ray flare was mea- 
sured with 10 (7 significance by EGRET at a flux level of 
(1.72±0.42) x 10" ^ cm" 2 s" 1 , 12 time s higher than the pre- 
vious upper limit (Bl oom et al.l Il997|) . The EC, in which 
the broad-line emission flux has been Comptonized, has been 
suggested for t he interpretation of the 7-ray emission i n the 
1997 flare (e.g. lMadeiski et al.lll999l:lRavasio et al.ll2002l) . 

In the VHE 7-ray range, the Crimean Ob servatory has 
claim ed a detection with 1.2 a significance (Neshpor et al. 
120011) . while HEGRA, observing in the same peri od, ob- 
tained only a much lower upper limit ( Kranich 2003, see de- 
tails in Discussion). Other p ast observations of this target 
resulted in upper limits only dC atanese et al. 1997]; iKranichl 
120031: lAharonian et al.ll2004l:lHoran et al.ll2004l) . 

In this Letter we report about the discovery of VHE 7-ray 
emission from BL Lacertae in 2005. Simultaneous observa- 
tion in the optical band in 2005 and 2006 permitted a search 
for correlations between the optical and VHE 7-ray activities. 

2. OBSERVATIONS AND DATA ANALYSIS 
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Fig. 2. — ALPHA-distribution of the 2005 data. The filled circles represent 
ON data. The light crosses correspond to normalized OFF data and a dotted 
curve describes a second order polynomial fit to the distribution of the OFF 
data. The vertical line indicates the ALPHA selection condition, which yields 
a total excess of 216 events at a significance level of 5.1 a. 

BL Lacertae (R.A.22 h 02 m 43.3 s , decl.+42°16'40"; J2000.0) 
was observed with the 17m diameter MAGIC (Major At- 
mospheric Gamma Imaging Cerenkov) telescope on the Ca- 
nary Island of La Palma (N28.2°, W17.8°.8, 2225 m a.s.L). 
The MAGIC telescope is the world's largest Imaging Atmo- 
spheric Cerenkov Telescope (I ACT). The telescope parame- 
ters and performance ar e described in iBaixeras et al.l (2004) 
and lCortina et al.l (120051) . The source was observed for 22.2 
hr from 2005 August until December. The telescope was 
pointing directly onto the object, recording so-called ON data. 
The background was estimated from observations of regions 
where no 7-rays are expected, which we define as OFF data, 
which were taken with sky conditions similar to those of ON 
data. 

The raw data were first calibrated ( Ga ug et aT]|2005h . and 
then processed wit h the standard MA GIC analysis and recon- 
struction software (Bre tz et al]|2005l) . Data runs with anoma- 
lous trigger rates due to bad observation conditions were re- 
jected from further analysis. The remaining ON data corre- 
sponded to 17.8 hr, while the OFF dat a correspond ed to 57.2 
hr. Image parameters of the raw data ( Hillas 1985) were cal- 
culated and compared for the ON and OFF data in order to 
check their consistency; excellent agreement was found. The 
selected data samples before the 7/hadron separation are com- 
pletely dominated by hadron events and the expected admix- 
ture of 7-ray events are below the statistical fluctuations of the 
data. 

In order to reject most hadronic background while still 
conserving the majority of the 7-candidates, a multi-tree 
classifier algorithm based on the "Rand om Forest (RF)" 
method (iBreimanl l200ll: iBock et al.ll2004l) was used for the 
7/hadron separation. The selection condition s were trained 
with Monte Carlo sim ulated 7-ray samples ( Kna pp & Heckl 
l2004t iMajumdar et al.l l2005|) and a sample of experimental 
background events. In the RF method for each event the so- 
called HADRONNESS parameter (H) was calculated from a 
combination of all the image parameters except the ALPHA 
parameter. ALPHA is the angle between the shower image 
principal axis and the line connecting the image center of 
gravity with the camera center. H assigns to each event a 
number between and 1 of being more hadron-like (high H) 
or 7-ray like (low H values). The selection of events with low 
H value enriches 7-ray events in the surviving data sample. 
Finally, the image parameter ALPHA was calculated and the 
ALPHA distribution for the ON and OFF data compared af- 
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Fig. 3. — Light curve of VHE 7-ray (> 200 GeV) flux as measured with 
the MAGIC telescope during 2005 [top]. A dotted horizontal line represents 
the average flux. Light curve of optical (R-bmd) flux during 2005 [middle] 
and 2006 [bottom] as measured with the 1.03m Tuoria and the 35 cm KVA 
telescope. In the optical light curve, the filled points represent the optical flux 
when simultaneous MAGIC observations were carried out. The measured 
average flux of those points is 9.2 mJy for 2005 and 4.2 mJy for 2006 



ter normalizing them on the area between 20° and 80° . Any 
7-ray signal should show up as an excess at small ALPHA 
values. The number of background events in the ON distri- 
bution at small ALPHA was determined by a second order 
polynomial fit (without linear term) to the ALPHA distribu- 
tion of the normalized OFF data. For our anal ysis we used an 
ALPHA cut of 9° (see also lAlbert et alll2006l l2007^bh . 

The cuts in H and ALPHA were optimized using data sam- 
ples from Crab Nebula observations at comparable zenith an- 
gles. The significance of an excess visible in the ALPHA 
distri bution was calculated according to Eq.17 in iLi & Mai 
(1983). The energy of the 7-ray events was also reconstructed 
by means of the RF method with MC 7- samples. The average 
energy resolution is estimated to be about 24% rms. 

Optical 7?-band observations were provided by the Tuoria 
Observatory Blazar Monitoring Program 23 with the 1.03 m 
telescope at the Tuoria Observatory, Finland, and the 35 cm 
KVA telescope on La Palma, Canary Islands. The magnitudes 
were then converted to linear fluxes using formula F[Jy] = 
3080.0 x i()-( ma g/ 2 - 5 ). 

Radio observations were also performed with UMRAO 24 at 
4.8, 8.0 and 14.5 GHz, and at 37 GHz with the Metsahovi 
Radio Observatory. 

3. RESULTS 

Figure Q] shows the local sky map for the 2005 7-r ay can- 
didates using the DISP method dLessard et al.ll200l|) . Only 
events > 350 photo-electrons were used, corresponding to an 
energy threshold of about 200 GeV 7-ray energy. The map 
was produced from the excess events distribution smoothed 
with a 2-D Gaussian of 0.1°. The black cross marks the nom- 
inal position of BL Lacertae. The small offset and the exten- 
sion of the image are comparable to the telescope point spread 
function (PSF, 0.1°) and the telescope pointing error (1.5'). 

23 more information at http://users.utu.fi/kani/lm/ 

24 UMRAO is partially supported by a series of grants from the NSF and 
by funds from the University of Michigan. 
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Fig. 4. — Differential energy spectrum of the 2005 BL Lacertae data. The 
solid line represents a power-law fit to the measured spectrum. The fit param- 
eters are listed in the figure. For comparison, the line fitted to the measured 
MAGIC CRAB spectru m using a power law with a changing photon index is 
shown by a dashed line (Wagner et al. 2005). 



In Figured the ALPHA distribution is shown. An excess of 
216 events over 1275.6 normalized background events yields 
a significance of 5.1 a for data above 350 photo-electrons. 
Two mostly independent analyses confirmed the result. 

The 2005 VHE 7-ray and the 2005 and 2006 optical light 
curves are shown in Figure [3] No significant evidence of flux 
variability in VHE 7-rays was found in the 2005 data. The de- 
rived average integral flux is F(E > 200 GeV) = (0.6 ± 0.2) x 
10 -11 cm -2 s -1 (x 2 /dof = 16.3/15), which corresponds to 
about 3% of the Crab Nebula flu x as measured by the MAGIC 
telescope (Wagner et al.l 120051) . In the optical light curve, 
the contribution from the host galaxy (1.38 mJy) was sub- 
tracted. The optical light curve shows a flare around the end of 
2005 October (MJD -53,670). Also radio light curves at 37 
GHz (A. Lahteemaki 2006, private communication) and 14.5 
GHz (UMRAO) suggest some minor flaring activity starting 
in 2005 November. 

The reconstructed differential energy spectrum (Fig. [4]) is 
well described by a simple power law: 



dN 
dE 



= (1.9±0.5)x 10~ n 



300 GeV 



-3.6±0.5 
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TeV s cm 2 



Because of the relatively steep slope, the systematic errors are 
estimated to be ~ 50% for the absolute flux level and 0.2 for 
the spectral index. 

Follow-up observations were carried out from 2006 July 
to Se ptember for 26.0 hr, using the so-called wobble 
mode (iDaum et al.ll997|) . where the object was observed at an 
0.4° offset from the camera center; 25.0 hr of the data passed 
all quality selection for the analysis. The DISP method was 
applied for the reconstruction of the shower direction with the 
final cut on the 2 parameter (the squared angular distance 
between the nominal source position and the reconstructed 7- 
ray direction). However, no significant excess could be found 
in the 2006 data. 

4. DISCUSSION 

The Whipple 10 m telescope observed BL Lacertae for 39. 1 
hr in 1995 and derived a flu x upper limit above 350 GeV at 
3.8% of Crab (iHoran et alJl20Q4h . HEGRA derived an up- 
per lim it above 1.1 TeV at 28 % of Crab with 26.7 hr obser- 
vation (lAharonian et al.ll2004l) . These upper l imits are con- 
sistent with our results. On the other hand, Neshpor et al. 
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Fig. 5. — SED of BL Lacertae. Black filled circles represent simultaneous 
2005 data of KVA and MAGIC as well as radio data taken by UMRAO and 
Metsahovi. Gray color points d escribe old me asurements (see detail in the 
inlay). The lines are taken from Ravasio et al. ( 2002). The solid line repre- 
sents a one-zone SSC model for the 1995 data, the dotted line is produced 
with SSC and EC components for the 1997 flare data. 



(1200 ll) claim VHE 7-ray detection using the GT 48 tele- 
scopes of the Crimean Astrophysical Observatory from data 
taken in summer 1998. The reported integral 7-ray flux is 
F(E > 1 TeV) = (2.1 ±0.4) x 10" 11 cm" 2 s"\ which is two 
orders of magnitude higher than the extrapolated value from 
this study. During the same period, 1998 July to August, no 
significant signal was found by HEGRA, and their reported 
flux upper limit in th e same energ y band is 7 times lower than 
the Crimean result (Kranich 2003). The VHE 7-ray detec- 
tion reported by Neshpor et al. cannot be explained without 
a remarkably huge and a very rapid flux variation offset by 
a few hours in consecutive nights from the HEGRA obser- 
vation. In the case of a leptonic origin of the 7-ray emis- 
sion, such a flare would possibly coincide with high activity in 
the optical as in the outburst of 1997 July, when the increase 
in flux was observed both in the optical and X -ray to 7-ray 
bands dBloom et al.lll997HTanihata et al.ll20QQh . However, no 
increased optical activity was detected during the Crimean ob- 
servation period (m p = 13.5-14.6, while 13.0- 14.6 in 2005). 
Throughout the EGRET observations for BL Lacertae, such 
a notably huge and rapid flare feature was never reported in 
the high energy 7-ray emission, which is considered to have 
the same origin as the VHE 7-ray emission regardless of the 
scenarios for the emission origin. 

In Figure [3j filled circles in the optical light curve repre- 
sent simultaneous observations with the MAGIC telescope, 



accepting ail day offset with respect to the MAGIC obser- 
vations. In 2005, there are 12 out of 16 nights with simultane- 
ous observations (average flux: 9.2 mJy), while 16 out of 23 
nights in 2006 (average flux: 4.2 mJy) have coinciding obser- 
vations. The absence of a significant excess of VHE 7-rays in 
the 2006 data indicates that the VHE 7-ray flux in 2006 was 
significantly lower than the flux level in 2005. In summary 
our results show similar tendencies both in the optical and the 
VHE 7-ray flux variations. Similarly, the 7-ray activity seen 
by EGRET observations in 1997 showed a strong correlation 
with optical activity. Such a correlation i s expected from lep- 
tonic origin scenarios (iBloom etal.119 97). Details of the 2006 
data and radio data will be discussed elsewhere. 

Figure \5\ shows the SED of BL Lacertae with results of 
this work and some previous data and model calculations by 
iRavasio et al.l (120021) . The VHE 7-ray points are corrected for 
the extra-galactic backgr ound light (EBL) ab sorption using 
the "Low" EBL model of iKneiske etafl (12004 . Our optical 
and VHE 7-ray points agree well with the solid line, which 
was derived using a one-zone SSC model on the 1995 data, 
whereas some deviations can be seen from the dotted line, 
which describes the 1997 flare data and involves SSC as well 
as EC components (IRavasio et al.ll2002h . To describe our re- 
sult such an additional EC component is not necessarily re- 
quired. 

BL Lacertae is the first LBL object with a clear detection 
of VHE 7-ray emission. The results of this work indicate that 
VHE 7-ray observations during times of higher optical states 
can be more efficient. Future long term monitoring of VHE 
7-ray emission could provide de tailed information o f a possi- 
ble periodicity predicted by e.g. lStirling et al.l (120031) and cor- 
relations with other wavelengths. Due to the observed steep 
spectrum, lowering the energy threshold of IACTs (e.g. the 
upcoming MAGIC-II project) would significantly increase the 
detection prospects for this new source class. 
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